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We analyze generalized Gaussian cat states obtained by superposing arbitrary Gaussian states, 
e.g., a coherent state and a squeezed state. The Wigner functions of such states exhibit the typical 
pair of Gaussian hills plus an interference term which presents a novel structure, as compared with 
the standard superposition of coherent states (degenerate case). We prove that, in any dimensions, 
the structure of the interference term is characterized by a particular quadratic form; in one degree 
of freedom the phase is hyperbolic. This phase-space structure survives the action of a thermal 
reservoir. We also discuss certain superpositions of mixed Gaussian states generated by conditional 
Gaussian operations or Kerr-type dynamics on thermal states. 
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I. INTRODUCTION 

The creation of quantum superpositions of living or- 
ganisms, what used to be merely a ridiculous possibility, 
may become feasible in a near future Of course, the 
center-of-mass states of a microbe, discussed by Romero- 
Isart et al jlj, are still very distant from the fantas- 
tic dead/alive Schrddinger cat @j. In any case, for the 
time being, one will have to be satisfied with inorganic 
quantum superpositions which, however, may involve a 
large number of particles ( "cat states" ) . At present they 
can be generated in the laboratory in a variety of sys- 
tems, e.g., optical cavities [3|, superconducting devices 
[U , Bose- Einstein condensates Q , free-propagating light 
beams [5], etc. Not only are such states crucial for testing 
the nonlocality of quantum mechanics outside the micro- 
scopic domain, but they also play an important role in 
some applications of quantum physics like information 
processing 0, 0] and quantum metrology 0-Q . 

The nonlocal quality of a cat state is most clearly dis- 
played when it is looked through the Wigner representa- 
tion. For instance, in the case of a superposition of two 
coherent states, the Wigner function exhibits two Gaus- 
sian "hills", corresponding to a classical superposition, 
and a nonclassical third Gaussian modulated by interfer- 
ence fringes halfway between the pair. This interference 
pattern is composed of alternating positive and negative 
straight bands parallel to the line joining the two Gaus- 
sian peaks [ToT | . 

The distributive property of the Wigner function al- 
lows one to construct the interference pattern of the su- 
perposition of several (or many) coherent states [TTj from 
the corresponding Wigner function for each pair. In fact, 
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as the set of coherent states is complete, any given pure 
state can be expressed as a coherent-state superposition: 
its Wigner function becomes a sum of cat-state Wigner 
functions. This observation is sometimes enough to un- 
derstand qualitative properties of general superpositions, 
like sensitivity to perturbations, or to the action of the 
environment (decoherence) [l2l. Il3j. 

In some proposed experiments the generated states 
consist of superpositions of coherent states lying on a 
manifold in phase space. This is the case, e.g., of the 
so-called arc states [l4], EH ■ The Wigner function of such 
a state is characterized by a banana-shaped positive re- 
gion decorated by interference fringes. Reciprocally, a 
state whose Wigner function is localized in the vicinity 
of a phase-space manifold, can be fitted by a combination 
of coherent states concentrated on that manifold [l6| . 

Consider now the evolution of a superposition of coher- 
ent states under a nonlinear unitary dynamics. In semi- 
classical regimes, for short times, an individual coherent 
state moves along a classical trajectory while suffering 
(approximately) a linear distortion. So, the wavepacket 
remains Gaussian but may become squeezed and rotated 
In general, different coherent states of the super- 
position will be distorted in different ways. This takes 
us to the following basic question: what happens to the 
Wigner function of a superposition of two coherent states 
when they are linearly, but independently, distorted? For 
instance, what is the interference structure of a superpo- 
sition of a coherent state and a squeezed state? 

This paper is dedicated to the study of the coherence 
structure of generalized Gaussian cats in the Wigner rep- 
resentation. In the pure case, a Gaussian cat is, by defini- 
tion, just a superposition of pure Gaussian states. Even 
though there is no natural way of defining what a su- 
perposition of mixed states would be, there are certain 
dynamics or protocols which produce what may be called 
mixed Gaussian cat states, given that the Wigner func- 
tion of such states are sums of Gaussians. We analyze two 
classes of mixed cats: (i) pure cats that suffered the de- 
cohering action of a thermal environment, and (ii) cats 
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generated either by conditional Gaussian operations or 
by a Kerr-like dynamics on input thermal states. In all 
cases we focus on the geometric structure of the interfer- 
ence fringes. 

The paper is organized as follows. We start by pre- 
senting our original dynamical motivation for studying 
generic superpositions of Gaussian states (Sect. [TTJ) . In 
Sect. IIIIl we summarize the geometrical formalism behind 
the Wigner representation which will be instrumental for 
deriving most of the results of this paper. In Sect. HVI we 
analyze in detail the problem of a superposition of arbi- 
trary pure Gaussian states. We show that for one-degree 
of freedom (two dimensional phase space) the interfer- 
ence pattern is in general hyperbolic (instead of the linear 
pattern exhibited by superposition of coherent states). 
For higher dimensions we derive a general normal form 
which shows that, in a suitable canonical coordinate sys- 
tem, the fringe pattern is still hyperbolic in each canon- 
ical plane. Curiously, this pattern survives the action of 
a linear Markovian environment (Sect. [V|, although its 
amplitude is progressively reduced over time. 

Section [VI] considers the result of using standard 
schemes for generating pure cat states, but, instead of 
injecting a coherent-state at the input, one feeds the sys- 
tem with a thermal state. This also produces mixed cat 
states which, however, are in general structurally differ- 
ent from the decohered Gaussian cats of Sect. IIVI 

Section [VIII presents our concluding remarks. 



II. 



A MOTIVATION FROM SEMICLASSICAL 
DYNAMICS 



A typical method for the approximate semiclassical 
propagation of quantum states consists of decomposing 
the initial state (possibly delocalized) into a suitable su- 
perposition of coherent states which are themselves prop- 
agated approxim ately by taking advantage of their qua- 
siclassical nature [18| . 

In the crudest scheme, each coherent state moves along 
a classical trajectory while preserving its shape (frozen 
Gaussian approximation). In the next level of improve- 
ment, an initial coherent state stays Gaussian during its 
evolution but may become squeezed and rotated (thawed 
Gaussian approximation) [19j . Both schemes correspond 
to approximating the exact Hamiltonian by its Tay- 
lor expansion around the instantaneous center of the 
wavepacket (q(t),p(t)) — x(t): 



H(x,t) = H(x(t)) + (x-x(t)) 
d 2 H 



dH 

dx 



+ (x-x(t)) 



dx 2 



(x- 



(t) 
x(t)) - 



(1) 



x{t) 



Here x represents the pair of canonical operators (q,f>). 
The frozen/thawed schemes are obtained by truncating 
the Taylor expansion to first / second order in x fPA fioj ] . 
In both cases the dynamics generated by H (x, t) is linear. 



Let us illustrate the semiclassical schemes above with 
a numerical example. Consider the propagation of a 
squeezed state in the Kicked Harmonic Oscillator (KHO). 
In suitable units the KHO Hamiltonian reads \2S : 



H{q,p,t)= l -(f + q 2 ) 



Kcos(q) 2^ 5(t-UT). (2) 



We choose an initial state \ipo) which is centered at the 
phase-space origin and squeezed along the g-axis. By ap- 
propriate choice of the parameters r and K, strong non- 
linear effects can be clearly seen after two pulses (kicks) . 
See Fig. [TJ 




FIG. 1: (color online) Nonlinear dynamics in the kicked har- 
monic oscillator. Shown are Wigner functions of (a) the 
evolved state at t = 2~; (b) the initial squeezed state (t = CP), 
which is centered at the origin but has ben displaced for clar- 
ity; (c) a coherent state, displayed for reference. Parameters 
are K = 2, r = tt/3, h = 0.0128. 

Our approximate description of the final state starts 
by decomposing the initial squeezed state \ipo) in terms 
of a one-parameter family of coherent states |</>(<?')) with 
centers lying on the line p = 0: 



dq'C(q')\^q')), 



(3) 



where (q\ipo), (q\(j)(q')) and C(g') are real Gaussians [21[ 
(in numerical calculations the integral is substituted by a 
finite sum). Then we propagate each \<p(q')} using the in- 
stantaneous local quadratic Hamiltonian, i.e., in Eq. ([2]) 
we make the approximation 



cos(g) fa a(t) + b(t) q + c(t) q 2 



(4) 



with a(t) = cos(q(t)), etc. In this way the final state 
becomes a superposition of Gaussian states, their centers 
threaded by a curved manifold in phase space (see Fig. [2]). 
Figure [3] presents a numerical comparison between ap- 
proximate and exact Wigner functions. The excellent 
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Fig. [T] can now be viewed as arising from the superpo- 
sition of the elementary patterns for each pair of final 
Gaussian states. In principle each wavepacket suffers a 
different distortion, so, the cat states obtained by super- 
posing pairs of wavepackets are in general nonstandard 
in the sense that their interference fringes are not par- 
allel. Inspection of Fig. ([2]) suggests that in the general 
case the interference fringes may be hyperbolic. 

In Sect. IIVI below we formulate and prove a precise 
statement about the hyperbolicity of the interference pat- 
tern of the superposition of two arbitrary pure Gaussian 
states. But first we need to review some mathematical 
results and notations which are extremely useful when 
dealing with Wigner functions and their related charac- 
teristic functions. 



FIG. 2: (color online) The initial squeezed state of Fig. [T] 
is decomposed into coherent states lying on the g-axis (top, 
schematic). Then each coherent state is evolved according 
to the local linear dynamics, so that it preserves its Gaus- 
sian shape but may rotate and get squeezed. Two examples 
are shown (a,b), together with their joint cat-state Wigner 
function. At t — 2~ the q-axis has evolved into the curved 
manifold. Same parameters as in Fig. [1] 



agreement indicates that the global nonlinear dynamics 
of the initial state can be understood, both qualitatively 
and quantitatively, as the collective effect of a swarm of 
localized states evolving linearly along classical trajecto- 
ries. 
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FIG. 3: Exact (thin line) vs. approximate (circles) Wigner 
functions at time t = 2~ . We exhibit the section q = —2 of the 
Wigner function, i.e., W(q = — 2,p) vs. p. Same parameters 
as in previous figures. 

The interference pattern of the Wigner function of 



III. MATHEMATICAL FORMALISM 

Most of the formulae of this section are widespread 
in the literature, sometimes under different notations 
[17L I'i'J 2 I . For the sake of self-completeness we present 
them compactly using a unified notation. The idea of the 
present formalism is to exploit the geometry of transla- 
tions and reflections behind the Wigner representation 
[lH . In this way, many of the cumbersome integrals usu- 
ally needed in Wigner calculus are reduced to simple al- 
gebra. 



A. Translations and Reflections 

Given a quantum state described by the density oper- 
ator p, its Wigner function is defined by 



W(x) 



1 



(2-Kh) n J_ 



+00 



dq' (q+\q'\p\q- 



-ip-q jh 



(5) 

We assume that the system possesses n degrees of free- 
dom. Canonical coordinates will be represented col- 
lectively by the column vectors q = (qi,...,q n ) T and 
p = (pi, . . . ,p n ) T (of course, T means "transposed"). A 
point in phase space is characterized by a 2rt-dimensional 
column vector x — (qi, . . . ,p n ) T . A similar notation is 
used for operators q, p, and x, e.g., q = (qi, . . . , q n ) T ■ (As 
an alternative to the mechanical point of view, one may 
adopt an optical perspective: If p represents the state 
of an n-mode quantum radiation field, then q and p are 
conjugate quadrature- vectors of the field.) 

Closely related to the Wigner function, the (symmetric 
ordered) characteristic function x(Q is defined as the ex- 
pectation value of the phase-space translation operator: 



with£= (Zj,tJ) T , and 



(6) 



(7) 
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In the last equation we have introduced the wedge (sym- 
plectic) product: 



(8) 



which can be turned into an ordinary scalar product (or 
matrix product) with the help of the symplectic matrix 
J, 



£ax = (J£) • x = i y x , 

where J is given by 

' On In 
In n 



(9) 



(10) 



Here l„ and 0„ are the n-dimensional identity and null 
matrix respectively. A symplectic matrix S is defined 
by the property of preserving the wedge product, what 
amounts to SJS T = J. 

The Weyl-Heisenberg translation operator is equiv- 
alent to Glauber's optical displacement operator D(a), 
which is expressed in terms of a complex parameter a 
and the annihilation/creation operators a/aJ, 



D{a) = exp(aa^ — a* a) . 



(11) 



The correspondence is established by making the iden- 
tifications a — (£ g + i^ p )/\2h and a = (q + ip)/V2h 

must 

Though not so widely known, the Wigner function can 
also be written as an expectation value: 



W{x) = 



1 



(irk) 7 



tr p R x 



(12) 



The (Grossman- Royer [24j |) operators R x are both uni- 
tary and hermitian. In particular, R x =o coincides with 
the usual parity operator, i.e., Ro\q) — \ — q). Given 
that the parity operator corresponds to a geometrical re- 
flection through the phase-space origin, and taking into 
account the property 



TxRqTI = R x , 



(13) 



the operator R x is interpreted naturally as the quantum 
version of the phase-space reflection through the point x 

Both sets of reflections and translations, {R x } and 
{Tj}, constitute orthogonal bases in operator space (with 
respect to the Hilbert-Schmidt product). Thus any op- 
erator A can be written as a linear combination of reflec- 
tions or translations 1231: 



A = 



1 



1 



dx A(x) R x , 



(14) 
(15) 



where the reflection and translation symbols, A(x) and 
respectively, are given by: 



A{x) = 2 n tr(i4 c ) , 
1(0 = tr(iff 



(16) 
(17) 



The symbols A(£) and A(x) are related via a (symplectic) 
Fourier transformation 



A(x) = 



1 



d£A(0e 



i £Ax/h 



(18) 



The relation between the Wigner and characteristic func- 
tions is a particular case ( A — p) of the general expression 
above: 

Wix) = j2^J d ^^ el ' AX/h - (19) 



B. Metaplectic transformations 

The quantum analog of a linear canonical (symplectic) 
transformation S is a metaplectic unitary operator Ms 
[I7I ]. We have already seen an example: the parity oper- 
ator Rq is the quantum version of the classical reflection 
S = — bn- A metaplectic operator can be thought as 
the quantum propagator associated with a Hamiltonian 
that is purely quadratic in the canonical operators q and 
P [13, E2I • The metaplectic operators together with the 
translation operators (associated with Hamiltonians that 
are purely linear in q and p) , constitute the set of unitary 
Gaussian operations, i.e., operations that leave the set of 
Gaussian states invariant [271 ] . 

Metaplectic operators respect the classical group com- 
position law, i.e., if S and S' are symplectic transforma- 
tions, then 



M S s' = ± M S M S > 



(20) 



where the ± s ign (unessential for our purposes) depends 
on both S,S' [l3,[2l- 

Metaplectics, translations and reflections interact ac- 
cording to the following formulas [13, HH : 



(21) 
(22) 
(23) 



These expressions entail the covariance of the Wigner 
function with respect to both translations and meta- 
plectic transformations. When a state is translated, 
p' = f^pT^, one deduces from Eq. (P2T1) that its Wigner 
function is also translated: 





- Rx+c 




= Tss, 


m s r x mI 


= Rsx- 



W p ,(x) = W p (x-0, 



(24) 
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using that — T_£, and the definition (|12p . Anal- 
ogously, for a metaplectically deformed state, p' = 
Ms/oMg, one easily proves the metaplectic covariance of 
the Wigner and the characteristic functions 



W p ,(x) 

xAZ) 



(25) 
(26) 



[using Ml = M s -i, and Eqs. (j6H2l22l23p ]. 

The reflection symbol Ms (a;) is an essential ingredient 
in our calculations of Sec. IIVI It is obtained from (fTO)) 
and given by [22H24j 



2 n i v fir- 
Ms(x) = — ; exp — x ■ C$x 



(27) 



where the symmetric matrix Cs is the Cayley transform 
of the symplectic matrix S 2^, 24j : 



C s = J (S - l an ) (S + l a „)" 



(28) 



An explicit expression for the integer v can be found in 



Ref. [24{. Like the ± sign in Eq. (|20p. v is a delicate 



object, and not really necessary for our dicussions. For 
these reasons we shall omit all details about this phase. 

Even when S + l 2n is singular, Eq. (f2"7| is still mean- 
ingful if we interpret that formula as a limit (in which 
one of the eigenvalues of S tends to —1) jl7j . 

In principle we could have worked in the q- 
representation. However, the expression for Ms(q, q') re- 
quires the explicit splitting of S into four blocks [ljj . This 
would lead to extremely cumbersome expressions when 
considering pairs of metaplectic transformations, like in 
Sec. HVl below. 



For the sake of compactness of forthcoming formulae, let 
us define a normalized Gaussian function: 



0(x;M,C) = ^|peKp[- (x-C)-M(x-Om , 

( 31 ) 

where x and £ are 2n-dimensional real vectors, and M is a 
symmetric 2nx2n complex matrix with positive real part, 
ReM > 0. The branch of square root is chosen in such a 
way that it reduces continuously to the positive root in 
the case of real M [28]. Note that if M is symplectic (a 
case we shall meet several times), then det M = 1 [29j . 

When Q (x; M, £) represents a probability distribution 
M _1 is the covariance matrix. Sometimes we shall abuse 
of the language and use such denomination for arbitrary 
Gaussians; in that cases "covariance matrix" is to be un- 
derstood as a shorthand for "the inverse of the matrix of 
the quadratic form 

The Wigner function of the fiducial state (|2"9"j) is readily 
calculated: 



W (x) = G (x;\,0) . 



(32) 



The Wigner function of the general Gaussian state |S, £) 
is obtained from the fiducial one by using the covariance 
properties (I24I25[) . The result is 



W lsx) (x)^g[x;(SS T y 1 X] ■ 



(33) 



In this case, the covariance matrix SS T is symplectic. 
It is worth showing also the corresponding characteristic 
function: 



X|s, O (e)=2-'^[e/2;(SS T )- 1 ,0] 



(34) 



Note that, for pure Gaussian states, W(x) coincides with 
x(2x), except for a factor and a shift of origin. 



C. Gaussian States 



IV. GAUSSIAN CAT STATES 



According to the standard definition, a state is said 
Gaussian if its Wigner function is a Gaussian. Any pure 
Gaussian state can be obtained from a fiducial Gaussian 
state by the combined action of a translation and a meta- 
plectic transformation [241 ] . 

For simplicity, we take as fiducial state the ground 
state of the n-dimensional isotropic harmonic oscillator. 
It will be denoted |0). In a system of units where fre- 
quency and mass are unity, uj = m = 1, the position 
wavefunction reads: 



(?|0) 



1 



(ttS)™/ 4 



-q 2 /2h 



(29) 



A general pure Gaussian state can be obtained from |0) 
by the successive application of a metaplectic operator 
and a translation: 



|S,C)=T C M S |0) 



(30) 



We shall call "Gaussian cat state" any superposition 
of two pure Gaussian states, i.e., 



I*) 



a\U,u) + b\V,v) , 



(35) 



where a and b are complex amplitudes, arbitrary to the 
extent that the state remains normalized. Its Wigner 
function reads 

W 9 (x) = \a\ 2 W lu . u) {x) + \b\ 2 W lv , v) {x) + \ab\2(x) . (36) 

The Wigner function of the superposition is the sum of 
the individual Wigner functions (two Gaussian "hills" 
centered at u and v) plus an interference term, given by 



X[x) 



(ttH) 71 



Re 



(U,«|^ \V,v) 



(37) 



with tp = arg(a*6). A typical example of the superposi- 
tion (|3"6")l was depicted in Fig. [5] 
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The purpose of this section is to explore the structure 
of the interference term. We start by writing the matrix 
element in the last equation as a vacuum expectation 
value: 



((J,u\R x \\/,v) = (0\M^R x f v M v \0) 



(38) 



First we use the fact that the composition of a translation 
and a reflection is also a reflection [22| . 



J£Ax/h ft 



+e/2 • 



(39) 
(40) 



in order to transform T^R X T V into a single reflection 
times a phase. After invoking Eqs. (|2QI23p we arrive at 

(U, U | J R x |V,«>=e faA ^+ i ^/ 2ft <0|A u - 1(:c _ 7)) M u - 1 v|0), 

(41) 

with 



C = u-v, 
rj = (u + v)/2. 



(42) 
(43) 



We sketch the next steps skipping the details: (i) Ex- 
pand My-iy into reflection operators with the help of 
Eqs. (|14l27p . (ii) Turn the resulting product of reflec- 
tions into a translation using the composition formula 

M 



R X Ry = e^^^ T 2 ( x -y) . 



(44) 



(hi) At this point one must calculate the average (0|Tg|0), 
for a certain £. But this is essentially the characteristic 
function of the vacuum, which can obtained from Eq. (|34p 
by setting S = I and £ = 0. (iv) Calculate the remaining 
Gaussian integral. The hnal result is 



|J2u-i(x-„)Mj-iv|0) = (TrK) n KG(x;G, V ) , (45) 



with 



K 



2 n i^ 



= , (46) 

^det[(U+V) + z(U-V) J] 

G = (UU T +VV T ) _1 [2-i(UU T - VV T )J] ,(47) 

and /i an integer. 

Summing up, the interference term may be written as 



X(x) = 2 \K\ Re 



ix AC / 'h-\-i<j> 



G{x;G,rj) 



(48) 



Here <j> stands for a phase that does not depend on x. Its 
precise value is not relevant for the forthcoming analysis. 

It can be checked that G is not only symmetric but is 
also complex symplectic, i.e., GJG T = J. Then, one has 
detG = 1 [H. 



A. Geometry of the Wigner Function 



The interference pattern (|48[) can be factored into a 
positive Gaussian envelope times an oscillatory function. 
The envelope is given by: 



-^env {x) 



2\K 



(nh) 



± exp [-(x -r,)-ReG(x- rj)/h] . (49) 



This Gaussian is centered at x = rj, the midpoint be- 
tween the centers of the individual Wigner distributions 
W\\j tU \(x) andW\v tV \(x). Its covariance matrix, (ReG) -1 , 
is an average of the individual covariance matrices: 



(ReG)" 1 = (UU T + VV T ) /2. 



(50) 



If U = V (equally distorted wavepackets), then the shape 
(covariance matrix) of the Gaussian envelope is equal to 
the individual Wigner functions. The height of the enve- 
lope depends on the amplitudes a and b of the Gaussian 
states [see Eq. (j35j) ]. e.g., for equal amplitude superpo- 
sitions, \a\ = \b\, the Gaussian envelope ([4*§f is twice as 
high as the individual Wigner functions. 
The oscillatory part reads 



Zosc(x) = cos [4> + x A Q/h + (x — rj) -IrnGfi-r)) /H] 



with 



ImG = - (UU T + VV T ) 1 (UU T -VV T )J. 



(51) 



(52) 



The first observation is that when U = V the quadratic 
term vanishes, and the phase becomes purely linear. In 
this well known case, the interference pattern is com- 
posed of straight lines parallel to the vector £ joining the 
centers u and v (see Fig.H]). 

Now we switch to the general case U ^ V. We want to 
characterize the quadratic form x ■ ImGx, which deter- 
mines the shape of the oscillations of the Wigner func- 
tion. This shape will be clearly identifiable when, by a 
suitable canonical coordinate transformation, we reduce 
ImG to a normal form. 

Without loss of generality, we shall assume V = bn. 
The general case can be recovered by setting U — > V _1 U, 
then making the canonical transformation x — > Vx, etc. 

Our starting point is the Euler decomposition [30] : any 
symplectic matrix may be put in a specific diagonal form 
by the use of two orthogonal symplectic transformations: 

S € Sp(2n, R) => S = OA s O' , (53) 

with 0, 0' G Sp(2n, E) n SO(2n) and 

A s = diag(si,...,s„,sj" 1 ,...,s^ 1 ) , > 1 . (54) 

We obtain as a corollary that SS T is diagonalized by a 
orthogonal symplectic matrix: 



SS T = OA ss rO T , 



(55) 
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FIG. 4: (color online) A gallery of Gaussian-cat-state Wigner functions. We show equal-weight superpositions, i.e., Eq. (|35p 
with \a\ = |6|. ^Frorn left to right: (a) Standard superposition of coherent states displaying a linear interference pattern, (b) 
Superposition of a squeezed state and a coherent state, (c) Superposition of two squeezed states both centered at the origin 
and with orthogonal squeezing directions. 



with the matrix of eigenvalues 

A ss t = diag(Ai,...A n ,A7/\...A~ 1 ) 



(56) 



where Aj = sf > 1. We are now in the position of calcu- 
lating the eigenvalues of ImG. With the choice V = \ 2n , 
Eq. ([52")) reduces to 

ImG' = - (UU T + l 2n ) _1 (UU T - l 2n ) J . (57) 
Using the corollary ([55]) as applied to UU T we get 

-0 T ImG' = AuuT ~ ' 2n 1 - ( E 0n \ >_fO n S 



Am it + I?,, V 0», 



l UU T T >2n 

where 5 is a diagonal matrix with positive entries: 

= - Xi ~ 1 = f) 
~* ~ Xi + 1 ~ 



(58) 



(59) 



The anti-diagonal matrix T ImG'0 can be diagonalized 
by the symplectic orthogonal matrix 



H= * ( '« !* 



(a 7r/4 rotation in each coordinate plane), that is, 
HO lmG'0H T = f „ S °'I 



(60) 



(61) 



So, we have shown that for any superposition (1351) there 
exists a canonical coordinate system (Q, P) where the 
quadratic part of the phase of oscillatory term ([BTjl re- 
duces to the normal form 



(62) 



with 0j > 0. In one degree of freedom, except for the 
degenerate case U = V, the pattern is always hyperbolic 
(see Fig. S]). 

The nonlinear interference patterns we described for 
continuous- variable systems may also be observed in dis- 
crete Wigner functions of finite-dimensional states. Con- 
sider, for instance, the phase-space representation of 
Grover's search algorithm for an A^-qubit system [31j . 
The computer starts in a pure momentum state (an equal 
superposition of all basis states) and, after some itera- 
tions, evolves into a position eigenstate corresponding to 
the searched item. At intermediate times the computer 
state is a weighted superposition of a momentum state 
and a position state. The Wigner function of such states 
(see Fig. 2 in Ref. [3l|) are very similar to the hyperbolic 
squeezed-state superposition in our Fig. [4] 32] . 



»=i 



V. DECOHERED GAUSSIAN CAT STATES 



This section discusses what happens to the Wigner 
function of a pure Gaussian cat state when it evolves 
under a general linear dynamics, i.e., a dynamics that 
preserves Gaussian states [13] ■ In the unitary case such 
dynamics are generated by quadratic Hamiltonians. We 
have seen that this corresponds to a metaplectic evolu- 
tion operator which transforms the Wigner function ac- 
cording to the covariance rule (|25p (if the Hamiltonian 
contains a term linear in x there will be some additional 
translation of the Wigner function). 

The most general linear evolution of a density operator 
p is described by the master equation (written in the 
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Lindblad form) |33{ : 



dp__ l_ 
dt ih 



H,p 



1 

2h 



M 

^2 (llL k p + pL\i k - 2L k k 



fc=l 



(63) 

with a quadratic Hamiltonian H and linear Lindblad op- 
erators Lfc, 



Lk = X k Ax 



H = ±x> Ex, 



(64) 
(65) 



Here B is a real symmetric matrix and a complex 
vector (for simplicity we assume that both H and L k are 
time independent). 

The corresponding evolution equation for the Wigner 
function W(x,t) is obtained by using the standard recipes 
for transforming master equations into partial differential 
equations [26ll34|. In the case of a quadratic (in x) Lind- 
blad equation we arrive at a so-called linear Fokker-Plank 
equation: 



dW 



d 



dx 



- (AxW) 



1 d 



2dx 



ox 



(66) 



The drift and diffusion matrices, A and D, respectively, 
can be compactly written as 



D = fiReT, 

A = J(B-ImT) , 



where 



M 



T 



A fc A 



(67) 
(68) 



(69) 



fc=i 



Note that D is real symmetric nonnegative, while, in gen- 
eral, A is just real (ImT is antisymmetric). In spite of 
some differences in the notation, the results above coin- 
cide with those in Ref. [33j (see also |35j). 

Let us now consider the evolution of the Wigner func- 
tion of an initially pure Gaussian cat state. The expres- 
sion for W(x, t) can be obtained by convolving the initial 
distribution W(xo, 0) [given by Eq. (|3"6"|) ] with the Fokker- 
Planck propagator: 



W(x,t) 



dx o P(x,t\x Oj 0) W(x o ,0) ■ (70) 



The propagator P(x, t\xo, 0), i.e., the Green function for 
Eq. is a real Gaussian function of both variables x 
and Xo [33, 36]. If W(xq,0) is a sum of Gaussians, then 
W(x, t) will also be. 

For an initial Wigner function like that in Eq. f|36[) , the 
integral above can be computed explicitly to get the de- 
sired solution at any given time t. However, we are only 
interested in the general structure of the Wigner function 
of the evolved Gaussian cat state, which is determined by 
the purely quadratic part of the evolving Gaussians. For 



this reason we shall focus on the evolution of the covari- 
ance matrices of the individual Gaussians that compose 
the cat state; there are four of them: two real ones and 
a complex-conjugate pair. 

Given a Wigner function satisfying the Fokker-Planck 
equation (|66l) , its associated covariance matrix obeys the 
following equation of motion (36j : 



^ = AC + CA T + D . 

dt 



(71) 



We remark that here and below, the "covariance" C 
stands for the inverse matrix of the quadratic form of any 
of the four Gaussians composing the cat-state Wigner 
function. At t = these matrices are symplectic but, 
when the dynamics (|71[) sets in, the symplectic symme- 
try is lost. 

The solution satisfying the initial condition C(i = 0) = 
C is given by [37j 

C(t) = e At C e ATt + /* dt'e^-^ D (*"*') . (72) 



This equation allows us to evolve one-by-one the covari- 
ance matrices of the individual Gaussian terms of the cat 
state p)]). 

It is convenient to split Eq. (|72p into real and imaginary 
parts: 



ReC(t) = e A *ReC e AT *+ /* dt'e A ^ D e^*-*'^) 

Jo 

ImC(t) = e A *ImC e ATt . (74) 

If ImCo is nonsingular, then, evidently, ImC(t) remains 
nonsingular for all times. Moreover, Eq. ([74"]) says that 
ImC(t) is congruent to ImCo, so, the full signature of 
ImC(t) remains constant in time (this is Sylvester's law 
of inertia [38|; "signature" indicates a triplet of integers: 
the number of positive, negative and zero- valued eigen- 
values). Similar congruence arguments, and the fact that 
ReCo is positive definite and D is nonnegative, imply that 
ReC(t) remains positive for all times. 

Now we have collected the basic ingredients to prove 
the main result of this section: the equation of mo- 
tion (|T2")) preserves the signatures of both Re[C _1 (t)] and 
ImfC- 1 ^)], i.e., 



RcfC- 1 ^)] - Re (C J ) 
MC- 1 ^)] ~ Im(C - 1 ) 



(75) 
(76) 



the symbol ~ meaning "congruent to" . 

This result leads to the following implications for the 
structure of the Wigner function that evolves from (|36l) . 
Obviously, the terms that evolve from Wi[j tU \(x) and 
Wiy jV \(x), remain positive Gaussians for all times, their 
covariances evolving according to Eq. ([731 [in these cases 
one has ImC(t) = 0, Vt]. Concerning the interference 
term I(x), the preservation of the signature of ImC -1 (t) 
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means that the oscillation pattern will still be repre- 
sented by the normal form (|62p in some coordinate sys- 
tem. However, differently from the pure case, the new 
coordinates will not be canonical in general. Of course, 
the oscillation pattern is multiplied by a Gaussian win- 
dow. 

The proof of (|75l76p for the interference term starts by 
expressing the inverse of C(f) in terms of its imaginary 
and real parts (omitting the time dependence): 



RefCT 1 ) = 



ReC + ImC (ReC)~ ImC , (77) 



ImfC- 1 ) = [imC + ReC (ImC) -1 ReC 



(78) 



Let us analyze the first line ((77)) : ReC > 0, then 
(ReC)" 1 > 0, then ImC (ReC)~ X ImC > (Sylvester's 
law). As the sum of positive matrices is positive, then 
Re (C _1 ) is positive for all times. 

In the second line ([78)) we use that, ReC being positive, 
it has a positive square root: ReC = M 2 . Then, defining 
ImC = M _1 (ImC)M _1 , which has the same signature as 
ImC, we get 



Im (C^ 1 ) = IVT 1 ImC + (ImC 



(79) 



Noting that ImC + (ImC) - ImC , we obtain 

Im (C" 1 ) - ImC - ImC - ImC ~ Im (C^ 1 ) . (80) 

The last congruence is a consequence of the initial cor- 
relation matrix Co being symmetric symplectic: Cq 1 = 
JC J T . □ 



VI. SUPERPOSITION OF MIXED GAUSSIAN 
STATES 



In the previous section we considered a pure cat state 
which was subjected to the decohering action of a linear 
reservoir. Now we study the opposite situation: a pure 
Gaussian state is first decohered and afterwards used as 
input for a cat-generating protocol. In what respects the 
Wigner-function interference pattern, we shall see that 
the latter case is more general, as elliptical structures 
may also appear. 

Consider, for instance, the cat-generating scheme of 

Fig. m SH|. 

If the initial Gaussian state p is pure, i.e., p — \ip)(?p\, 
then the final state is also pure, p' = \ip')(ip'\, with 



(81) 



where the sign ± depends on the result of the qubit 
measurement, and 7V± is a state-dependent normaliza- 
tion factor. When U is a linear operation (a combination 
of a translation and a metaplectic unitary), the output 
state (|5T)) becomes a pure Gaussian cat state, exactly like 





Po 



H 



H 



U 



P 



FIG. 5: Circuit for generating Gaussian superpositions. The 
top line carries a single qubit initialized in the state |0). The 
bottom line is fed with a continuous variable Gaussian state 
po- The circuit is composed of two Hadamard gates (H) and a 
controlled-!/ linear operation. Measurement of the outgoing 
qubit produces a Gaussian superposition p at the bottom 
output. 



those described in Sec. IIVI However, if the initial state 
is Gaussian mixed, then the final state, 



p'=M ± (l±U)p (l±&), 



(82) 



may be called a mixed Gaussian cat [401 ] ■ A possible im- 
plementation of the circuit above for generating mixed 
superpositions uses two free-propagating quantum opti- 
cal modes interacting via a cross-Kerr nonlinear crystal 
40]. 

A formally similar result is produced when a single 
mode propagates in a medium exhibiting a Kerr nonlin- 
earity. In the interaction representation, the evolution is 
governed by the Hamiltonian [4l| 



H = jrY , 



(83) 



where n = cva is the number operator, and 7 an energy 
scale. At the revival time 



T = 



2irh 

7 



(84) 



an initial coherent state |ao) is perfectly reconstructed. 
For times equal to (p/v)T, with pjv an irreducible frac- 
tion, the evolved state consists of a superposition of v co- 
herent states lying on a circle of radius |ao| (fractional re- 
vivals) [42j]. However, this phenomenon is not restricted 
to coherent states, because at fractional-revival times the 
Kerr propagator becomes a sum of harmonic oscillator 
propagators, i.e., 



-27ri fin 2 /v 



c fc M k , (85) 



where c k are complex Fourier coefficients having equal 
moduli [43| . The notation M k emphasizes the metaplec- 
tic nature of the harmonic oscillator propagators, which 
correspond to phase-space rotations of angles 2irk/v. 
Equation ()85|) shows that at revival times any state will 
be transformed into a superposition of v (rotated) repli- 
cas of itself. If the initial state po is Gaussian mixed, say 
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a displaced thermal state, it will evolve into the mixed 
Gaussian superposition 



E 

k,j=i 



c k c* M k po M) . 



(86) 



Figurc^a) exhibits an example for v = 4 (a mixed "com- 
pass" [12| state). Two kind of interference patterns can 
be identified. The patterns corresponding to opposite 
replicas are linear, while those corresponding to contigu- 
ous replicas are elliptical (circular). In all cases, the ex- 
tension of the interference regions are smaller than those 
corresponding to the diagonal terms M k Po M^. Had we 




p - 





-5 - 



-5 





q 



FIG. 6: (color online) Wigner function of a mixed compass 
state. An initial displaced thermal state is evolved with the 
Kerr Hamiltonian up to a fractional revival time (y = 4). 

started with a suitable thermal squeezed state we should 
have obtained hyperbolic fringes instead of elliptical ones 
(graphics not shown). 

Note the formal similarity between both cat-generating 
schemes described above: The interference terms gener- 
ated with the Kerr Hamiltonian [Eq. ((86)) ] can also be 
obtained from the scheme of Fig. [5] [Eq. (|82|)]. by us- 
ing alternatively 7r/4 or 7r/2 controlled rotations. So, in 
general, both schemes produce similar mixed cat states. 
Such cat states arc structurally different from the deco- 
hered cats of Sect. IVI1 which can only exhibit hyperbolic 
fringes in their Wigner functions. 

We shall not present a detailed analytical description 
of the Wigner functions of the mixed cat states (|82I86I) . 
Explicit expressions can be worked out along the lines of 
Sect. HV] (see also Ref. j45|). However, we would like to 
exhibit a simple analytical explanation of the shrinking 
of interference patterns (see Fig. [5]). 

Consider for simplicity the binary cat produced by the 



Kerr dynamics at half the revival time (y = 2): 

p' = \{l + iMn) Po (1 - iMl) = ~(1 + iRo) po (1 - iR ) ■ 

(87) 

Here M v denotes the half-period harmonic evolution, 
which is equivalent to the parity operation Rq, i.e., the 
reflection through the phase-space origin. The Wigner 
function is calculated from Eq. (TT21) : 

27rfi W'{x) = tip R x + trRopoRoRx + 2RcitrR p Rx ■ 



Using the cyclic property of the trace, and the composi- 
tion formulae for reflections and translations (|39I40I44[) . 
we obtain the cat Wigner function, 

2W'{x) = W (x) + W (-x)+4Reixo(-2x), (89) 

in terms of the Wigner and characteristic functions of the 
initial state po [44]. If po is a displaced thermal state, 
Po = T v pthT^, then we arrive at 

2W'(x) = W & (x+7i) + W t h(x- rj) - 4 sin(xA 2r)) Xth(2x) . 

(90) 

Remarkably the interference pattern of the cat Wigner 
function is described by the characteristic function of the 
thermal state. When temperature grows Wth.(x) becomes 
wider. Then, the characteristic function must shrink, be- 
cause Wigner and characteristic functions are related by 
a Fourier transform. So, by increasing the temperature 
one can make the interference pattern as small as desired 
40]. 



VII. CONCLUDING REMARKS 

We have studied the Wigner functions of general su- 
perpositions of Gaussian states. For the pure case, we 
showed that the structure of the interference pattern is 
hyperbolic in general. This structure is robust against 
the action of a linear environment. We also analyzed two 
families of mixed Gaussian cat states which may also ex- 
hibit elliptic fringes. 

Our approach was geometric, qualitative. Anyway, the 
analytical tools presented here can in principle be used 
for quantitative purposes, like descri bing ho w cat states 
lose coherence in linear environments 45l-l48j , or to deter- 
mine the order of nonclassicality [4!| of the generalized 
cats. 



As far as we know, generalized cat states, i.e., show- 
ing nonlinear intereference patterns in their Wigner func- 
tions, have not yet been created in the laboratory. Some 
squeezed superpositions have already been produced [f| , 
others may become reality soon (50l. [5lj[. However, such 
squeezed cats are still degenerate: as there is no relative 
squeezing between the superposed states, the Wigner in- 
terference pattern is linear. 

We would like to conclude by mentioning some theoret- 
ical generation of "hyperbolic" cats [32j . A set of quan- 
tum states localized on the classical periodic orbits of a 
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chaotic map, can be used as a basis in which the descrip- 
tion of the eigenstates of its quantum version is greatly 
simplified. This set can be improved with the inclusion 
of short time propagation along the stable and unsta- 
ble manifolds of the periodic orbits [52]. These "scar 
functions" [53[ , when viewed through a phase space rep- 
resentation look very much like the hyperbolic cat in our 

Fig. He) m. 
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